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(Dated: February 3, 2008) 

We report on precision measurements of the elastic cross section for electron-proton scattering 
performed in Hall C at Jefferson Lab. The measurements were made at 28 distinct kinematic 
settings covering a range in momentum transfer of 0.4 < Q 2 < 5.5 (GeV/c) 2 . These measurements 
represent a significant contribution to the world's cross section data set in the Q 2 range where 
a large discrepancy currently exists between the ratio of electric to magnetic proton form factors 
extracted from previous cross section measurements and that recently measured via polarization 
transfer in Hall A at Jefferson Lab. This data set shows good agreement with previous cross section 
measurements, indicating that if a here-to-fore unknown systematic error does exist in the cross 
section measurements then it is intrinsic to all such measurements. 



I. INTRODUCTION momentum transfer, Q 2 > 1 (GeV/c) 2 . This is due pri- 

marily to recent measurements from Hall A at Jefferson 

Recently, there has been much renewed interest in the 
proton electromagnetic form factors in the region of four- 
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FIG. 1: (Color online) Single photon exchange diagram for 
electron proton elastic scattering. 



Lab [1|, |2| on the ratio of the Sachs electric to magnetic 
form factors via the polarization transfer technique 0,13 • 
These data are in stark disagreement with previous ex- 
tractions of these form factors 0, B @ from cross sec- 
tion measurements utilizing the Rosenbluth separation 
technique jlOj . 

There have been recent efforts ^2 E2- El to extract 
the individual form factors by combining the cross section 
and polarization transfer results. However, it is clear that 
the data sets from these two techniques are systemati- 
cally inconsistent ^2] an d, as such, the method that one 
chooses for combining the data sets is not well defined. 
It is critical, then, that the source of the discrepancy be 
identified, if there is to be any chance of pinning down 
the Q 2 dependence of the individual form factors |l3( . 

In this paper we will present results from 28 new pre- 
cision measurements of the ep elastic cross section in the 
range 0.4 < Q 2 < 5.5 (GeV/c) 2 performed in Hall C at 
Jefferson Lab. Although the kinematics are such that 
only limited Rosenbluth separations of the form factors 
can be performed, this data represents a significant con- 
tribution to the world's cross section data set, and as 
such, can help provide new constraints on global fits from 
which the form factors can be extracted. 

The high precision and large kinematic coverage of the 
new Jefferson Lab data can help provide crucial informa- 
tion as to whether there exists an experimental system- 
atic error in the world's cross section data set, which is 
dominated by the data from SLAC in the Q 2 range where 
the discrepancy with the polarization transfer data exists. 



II. EP ELASTIC SCATTERING 

The elastic scattering of an electron from a proton tar- 
get can be represented in the first order Born approxima- 
tion by the exchange of a single virtual photon between 
the leptonic and hadronic electromagnetic currents. This 
exchange is represented by the diagram in Fig. ^ and is 
often referred to as the one photon exchange approxima- 
tion (OPEA), with 4- momentum transfer 



k'„ 



(1) 



where k p (k') is the four-momentum of the electron be- 
fore (after) scattering. For space- like photons (q 2 = q^q^ 
< 0) it is customary to define the absolute value of the 
square of the four-momentum transfer 



Q 2 = -q 2 « 4EE' sin 2 (9/ 2). 



(2) 



If the proton were point-like, then the cross section 
could be calculated within the framework of quantum 
electrodynamics (QED) to give 



_ da pl _ E' a 2 cos 2 (9/2) 



dfl E 4£ 2 sin 4 (6>/2)' 



(3) 



However, the spatial extent of the electromagnetic charge 
and current densities of the proton lead to the introduc- 
tion of form factors, which modify the proton vertex and 
parameterize the protons internal structure. It is com- 
mon to see the cross section expressed in terms of the 
Sachs electric and magnetic form factors, Ge p and Gm p - 
These form factors are defined in such a way that only 
terms quadratic in them appear in the Rosenbluth ex- 
pression for the cross section, 



da 

dh 

where 



G 2 Ep (Q 2 )+rGj Ip (Q 2 ) 

1+T 



2rG 2 Mp (Q 2 )tan 2 (9/2) 



(4) 
In 



r = Q 2 /AM 2 and M p is the proton mass 
the non-relativistic limit, Ge p is given by the Fourier 
transform of the spatial charge distribution, while Gm p 
is given by the Fourier transform of the spatial magne- 
tization distribution. At zero momentum transfer, the 
proton is resolved as a point particle of total charge equal 
to one and total magnetic moment given by [x p = I + n p , 
where k p — 1.7928 is the proton anomalous magnetic 
moment. This leads to the normalizations 



Ge p (0) = 1 and G Mp (0) = /V 



(5) 



III. 



EXTRACTION OF FORM FACTORS FROM 
CROSS SECTION MEASUREMENTS 



The Rosenbluth expression, Equation can be recast 
in terms of the relative longitudinal polarization of the 
virtual photon, e = [l + 2(1 + r) tan 2 (6>/2)] -1 , as 



da 

dh 



eG 2 Ep (Q 2 )+rG 2 Mp (Q 2 ) 



e(l + r) 

with the reduced cross section defined by 
_ da s(l + t) 



dn 



a pi 



eG 2 Ep (Q 2 ) + rGi Ip (Q 2 



(6) 



(7) 



At fixed Q 2 , the individual form factors, Ge and Gm, 
can be extracted from a linear fit in e to the measured 
reduced cross section. Such a fit is generally referred to 
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FIG. 3: Schematic of the Hall C beamline. 



backgrounds were found to be less than 0.1%. 



FIG. 2: Fractional contribution of Ge to the cross section 
assuming Gm/^p = Ge (form factor scaling). 



Hall C Beamline 



as a Rosenbluth fit and yields tG 2 m as the intercept and 
G\ as the slope. Due to the r weighting of G 2 M , the 
cross section becomes less sensitive to Ge at large Q 2 . 
Hence, the accuracy with which Ge can be extracted de- 
creases (inversely) with Q 2 and Rosenbluth separations 
eventually fail to provide information on the value of Ge- 
This failure was part of the impetus for the develop- 
ment of the polarization transfer technique. The frac- 
tional contribution of Ge to the cross section assuming 
Gm/ Hp = Ge (form factor scaling) is shown as a func- 
tion of e in Fig.Elfor Q 2 values of 1, 3, and 5 (GeV/c) 2 . 
At Q 2 = 3 (GeV/c) 2 , G E contributes only 12% to the 
cross section at e — 1, with this contribution decreasing 
approximately linearly as e — > in this Q 2 range. 



IV. EXPERIMENT 

The ep elastic scattering data presented here were ob- 
tained as part of experiment E94-110 |2^|. which was 
intended to separate the longitudinal and transverse un- 
polarized proton structure functions in the nucleon res- 
onance region via Rosenbluth separations. The exper- 
iment utilized the high luminosity electron beam pro- 
vided by the CEBAF accelerator and was performed in 
Jefferson Lab Hall C during summer and fall of 1999. 
Scattered electrons were detected in the High Momen- 
tum Spectrometer (HMS). Additionally, the Short Orbit 
Spectrometer (SOS) was used to detect positrons, which 
were used to determine possible electron backgrounds 
originating from charge-symmetric processes such as tt° 
production and subsequent decay in the target. For the 
kinematic of the elastic scattering measurements, these 



The Hall C beamline from the beam switch yard to 
the beam dump in the experimental area is shown in 
Fig- El The beam from the accelerator south linac enters 
the Hall C Arc and passes through a series of dipole and 
quadrupole magnets which steer it into the Hall. The 
beam position and profile can be measured at several 
stages in the Arc with the use of superharps. The super- 
harps consist of a set of fine wires (two horizontal and 
one vertical) which are moved back and forth through 
the beam to determine the centroid position to about 
10 /im. However, these measurements are invasive and 
can not be performed during data taking. Continuous 
monitoring of the beam position in the Arc is done with 
the aid of three beam position monitors (BPMs), which 
are nondestructive to the beam and are calibrated with 
superharp scans. 

The absolute beam position provided by scans of each 
of the three superharps allows the trajectory of the beam 
through the magnets to be determined. This, combined 
with knowledge of the field integrals of the Arc mag- 
nets, then allows the absolute beam energy to be deter- 
mined to better than 0.1%. Absolute beam energy mea- 
surements which require superharp scans were performed 
about twice per beam energy setting. 

Accelerator cavity RF instabilities have been observed 
to cause variations in the beam energy of about 0.05%. 
These variations of the beam energy can be measured 
using the relative positions provided by the Arc BPMs. 
These BPMs were read into the data stream every second 
and used to monitor the beam energy drift. In principle, 
the effect of a drift can be corrected for if a large enough 
sample of events is considered. However, the effect of 
beam energy drift on runs used in the current analysis 
was studied and found to be less than a 0.02% effect on 
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the beam energy. 



The beam position monitoring system in the Hall con- 
sists of three BPMs and two superharps for calibrations. 
Deviations in the angles of the beam on target translate 
into corresponding offsets in the reconstructed angles, 
whereas deviations in the vertical (spectrometer disper- 
sive direction) position of the beam will manifest them- 
selves a s apparent momentum and out-of-plane angle 
offsets in the spectrometers. The effect of a beam po- 
sition offset can be calculated from the optical matrix 
elements for the spectrometer. For a 1 mm vertical off- 
set of the beam on target, the shifts in the reconstructed 
momentum and out-of-plane angle in the HMS are about 
0.08% and 1 mrad, respectively. 

The centroid of the beam spot, determined by the 
beam steering into the Hall, is constantly monitored by 
both fast-feedback electronics and visual displays of the 
BPM readouts and is adjusted to prevent large drifts of 
the on-target position during data taking. A study of 
the run-to-run beam steering stability was made during 
the running of this experiment. In this study, the run- 
to-run variations in the vertical position on target were 
measured to be less than 0.2 mm, resulting in a corre- 
sponding point-to-point uncertainty in the reconstructed 
momentum of 0.0f6%. The run-to-run variations in the 
angles on target were found to be less than 0.04 mrad. 

In order to minimize localized target boiling effects in 
the liquid hydrogen, the small intrinsic beam spot size of 
about 300 /im 1 was increased by a set of fast rastering 
magnets before entering the Hall. The fast raster pro- 
duced a rectangular pattern with a full width of about 
4 mm in the horizontal and 2 mm in the vertical. Cor- 
rections due to the vertical rastering were calculated and 
corrected event- by- event. 

The beam current monitoring system in the Hall con- 
sists of two (three for E94-1 10) resonant microwave cavity 
beam current monitors (BCMs) . The BCMs provide con- 
tinuous measurement of the current and are calibrated 
to about 0.2 /i A by use of an Unser monitor in the Hall. 
Dedicated calibration runs were performed about once 
every three days during this experiment to minimize the 
effects of drifts in the BCM gains. The current was care- 
fully monitored during data taking and was required to 
be 60 ± 2/iA. The normalization uncertainty due mostly 
to the Unser was estimated to be 0.4%. The run-to-run 
uncertainty in the beam current of 0.2% was estimated 
by combining in quadrature the fit residuals from the 
calibration runs and the typical observed drift between 
calibrations. Detailed information on the current moni- 
toring systems in Hall C can be found in Reference fisj ]. 
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[1] This is about 3x the normal intrinsic spot size. 



FIG. 4: Representation of the cryogenic and dummy target 
assembly. 



B. Target 

A representation of the cryogenic target assembly is 
displayed in Fig. 0] and shows the three "tuna can" 
shaped cryogen cells, as well as the dummy target. Each 
can was machined out of aluminum to provide a very uni- 
form cylindrical shape which 'bulges' a negligible amount 
when the cell is pressurized to about 25 psia [l9|. The 
hydrogen cell was measured to have an inside diameter of 
40.113 mm when warm and 39.932 mm when cold, and 
a cylindrical wall thickness of 0.125 mm. Due to the cir- 
cular shape, the average target length seen by the beam 
depended upon both the central position of the beam 
spot and the size and form of the raster pattern. The 
normalization uncertainty in the hydrogen target length 
was estimated to be 0.3% and the run-to-run uncertainty 
was estimated to be 0.1%. The dummy target was made 
from two 0.975 mm thick, rectangular sheets of aluminum 
separated by 40 mm. Additional details on the target as- 
sembly can be found in Reference |l9j . 

Localized target density fluctuations, induced by an in- 
tense incident beam, can modify significantly the average 
density of a cryogenic target. Uncertainties in target den- 
sity enter directly as uncertainties in the total cross sec- 
tion, and can be current-dependent on a point-to-point 
basis. The current-dependence can be measured by com- 
paring the yields at fixed kinematics with varying beam 
currents. The deadtime-corrected yields should be pro- 
portional to the luminosity (and, therefore, to the target 
density). 

The result of such a 'luminosity scan' for E94-110 is 
shown in Fig. [5] where the luminosity relative to the 
lowest current has been plotted on the vertical axis. The 
error bars on the data are statistical only and do not 
reflect fluctuations in the beam current. The correction 
factor applied to the measured target density (at zero 
current) to account for the reduction resulting from lo- 
calized target boiling is given by the product of the fitted 
slope and the current at which the data was taken. For 
the present data the current was typically kept at 60 ± 
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FIG. 5: Relative hydrogen target yield versus beam current. 

2 [iA, resulting in a density correction of 

(2.4 ± 0.2)% ■ 60/iA/100/iA = 1.44 ± 0.12%. (8) 

The uncertainty in the current did not contribute appre- 
ciably to the uncertainty on this correction. The total 
estimated run-to-run uncertainty in the target density is 
0.1%. 

C. HMS Spectrometer 

The HMS is a magnetic spectrometer consisting of a 
25° vertical bend dipolc magnet (D) for momentum dis- 
persion and three quadrupole magnets (Ql, Q2, Q3) for 
focusing. All magnets are superconducting and were op- 
erated in a mode to provide a point-to-point optical tune. 
A schematic side view of the HMS is shown in Fig. El and 
includes representations of the pivot (with target cham- 
ber), magnets, and the shielded hut containing the de- 
tector stack. 

The detector stack is shown in Fig. [7| and consists of 
two vertical drift chambers [2(j (DC1 and DC2) for track 
reconstruction, scintillator arrays (S1X(Y) and S2X(Y)) 
for triggering, and a threshold gas Cerenkov and elec- 
tromagnetic calorimeter, which were both used in the 
present experiment for particle identification (PID) and 
pion rejection. 

The acceptance limits of the HMS in-plane (Y 1 ) 
and out-of- plane (X 1 ) scattering angles are defined by 
an octagonal collimator positioned between the tar- 
get and the first quadrupole magnet. The edges of 
this collimator define a maximum angular acceptance of 
-28 < Y' < 28 mrad and -75 < X 1 < 75 mrad, and 
a total solid angle of about 6.8 msr. Additional details 
on the HMS can be found elsewhere 261. 



FIG. 6: (Color online) Schematic drawing of the HMS spec- 
trometer. 
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FIG. 7: Schematic drawing of the HMS detector stack. The 
component detectors are described in the text. 



D. Data Acquisition 

Data acquisition was performed using the CEBAF On- 
line Data Acquisition (CODA) software [24[ running on a 
SUN Ultra-2 workstation. The detector information for 
each event was collected from the front-end electronics 
by VME/CAMAC computers (collectively referred to as 
Read-Out Controllers or ROCs). Event fragments from 
the ROCs were then transfered via TCP/IP to CODA, 
which formed events and wrote them to disk. 



V. DATA ANALYSIS 

CODA events from individual run files where decoded 
by the Hall C Replay software, which reconstructed the 
trajectories of individual particles from hit information 
in the drift chambers. Tracks were then transported back 
to the target via an optical transport model of the HMS, 
which allowed the determination of the particle kinemat- 
ics. For each run an HBOOK [2^| ntuple was then created 
which contained the reconstructed event kinematics and 
calibrated PID detector information. The final analysis 
of the ntuples into experimental yields is described in the 
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sections which follow. 



A. Kinematic Calibrations 

One of the larger e-dependent uncertainties that di- 
rectly affects Rosenbluth separations is that due to the 
uncertainties in the kinematics at which the cross sections 
are measured. It is convenient to absorb this uncertainty 
directly into the cross sections by calculating the ex- 
pected difference in the measured cross section when the 
kinematics are changed from the nominal values within 
their uncertainties. In order to minimize this uncertainty, 
it was critical that the kinematic quantities, E, E', and 
9 be determined to the best possible precision. This was 
aided by the kinematic constraint of elastic scattering, 
that the reconstructed mass of the unmeasured hadronic 
state be equal to the proton mass. 

For each kinematic setting, the difference of the recon- 
structed invariant mass, W, from the proton mass (AW 
— W - M p ) was calculated after correcting for the effects 
of energy loss due to both ionization and bremsstrahlung 
emission. This provided a large set of kinematics for 
which the dependence of AW on possible energy and an- 
gle offsets could be studied. Finally a minimization of 
AW was performed to determine the best set of kine- 
matic offsets under the following assumptions: 1) the 
offset of the nominal HMS central momentum from the 
true value was a constant fractional amount, and 2) the 
offset of the nominal HMS central angle from the true 
value was a constant. The nominal HMS momentum 
used in this study was that determined in 26], while 
the nominal HMS central angle was determined from a 
comparison of marks scribed on the floor of the Hall to a 
marker on the back of the spectrometer, which indicated 
the optical axis. 

The reconstructed W values for these data are plotted 
versus scattering energy in Fig.[S| for eight different beam 
energies and thirty one unique kinematic settings. It was 
found that the entire data set could be well described by 
assuming that the true HMS central angle was smaller 
than the nominal value by 0.6 mrad, and that the true 
HMS central energy was smaller than the nominal value 
by 0.39%. 

The true beam energy was also found to be smaller 
than the Arc measurements by an amount that varied 
with the energy. This latter result was subsequently con- 
firmed @| by a remapping and analysis of the field for 
one of the Arc magnets. The reconstructed values for W 
are shown, both before (open symbols) and after (solid 
symbols) correcting for the kinematic offsets found from 
these studies. The corrected values are all seen to be 
within 1-2 MeV of the proton mass. This procedure was 
used for the E94-110 data and yielded a estimated un- 
certainty in the corrected beam energy of 0.056%, about 
half that typically quoted from Hall C Arc measurements. 
We estimate that this uncertainty is the quadrature sum 
of equal normalization and run-to-run uncertainties. We 
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FIG. 8: (Color online) Reconstructed W vs. HMS central 
angle for elastic scattering kinematics. Open symbols repre- 
sent data before kinematic corrections were applied, while the 
solid symbols represent the data after applying the calibration 
corrections. 



note that the beam energies determined from the Arc 
measurements utilizing the updated field maps agree with 
the current results to about 0.05%. 



B. Binning the Data 



The data were binned on a 2-dimensional grid in the 
reconstructed variables E' and 0. This was because, at 
fixed beam energy, the inclusive cross section only de- 
pends on the scattered electron energy and angle. In 
practice, the binning in E' was converted to a binning in 
SP/P, which is the more natural variable for the applica- 
tion of the acceptance corrections. The ranges were cho- 
sen such that the entire angular acceptance was included 
and the 8P/ P acceptance was well determined from the 
model of the HMS. For 8P/P, the binning chosen was 16 
bins over a range of ±8%, while for 9 the binning chosen 
was 20 bins over a range of ±35 mrad (AO = 3.5 mrad). 
We note that the physical solid angle coverage (Ail) can 
be different for each AO. 
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C. Analysis Procedure 

For a beam of electrons of energy E incident on a fixed 
proton target, the number of electrons scattered at an 
angle 9 in a solid angle Ail is related to the differential 
cross section, da(9)/dVL by 



O) 



where C is the integrated luminosity. This is not the 
OPEA cross section of Equation but contains contri- 
butions from higher order QED effects. These include 
virtual particle loops, multi-photon exchange, as well as 
the emission of bremsstrahlung photons, both before and 
after the scattering. 

The emission of unmeasured bremsstrahlung photons 
by either the electron beam or the outgoing detected elec- 
tron results in energies at the scattering vertex which are 
either smaller (the incoming case) or larger (the outgo- 
ing case) from those used in the reconstruction of the 
kinematics. This results in a large radiative "tail" in 
both the reconstructed E' and the invariant hadron en- 
ergy distributions for elastic events. To compare to the 
OPEA cross section, requires that this radiative tail be 
integrated to some cut-off in E', with a correction fac- 
tor, which included the remaining higher order effects, 
depending on this cut-off. This "radiative" correction is 
applied as a multiplicative factor (denoted RC) and is 
discussed in more detail in section IVHl Because the ra- 
diative tail extends beyond the threshold for single pion 
production at W 2 « 1.16 GeV 2 , the integration was cut- 
off at W^ax < 1-16 GeV 2 to avoid including events from 
inelastic processes. The corresponding correction factor, 
RC(W^ ax ), is therefore cutoff dependent. 

In addition, the measured number of counts must also 
be corrected for detector efficiencies, eff, and the effec- 
tive solid angle acceptance, AD, e ff(8,E'), after subtrac- 
tion of counts from background processes, BG(0, E'). In 
this experiment the measured cross section was deter- 
mined for each bin on a 2-dimensional grid of the elec- 
tron scattering energy and angle, E' and 9, across the 
entire phase space for which the spectrometer has a non- 
zero acceptance. The extracted cross section was then 
determined from the relation, 



jg^yCg) = RC{Wl ax ) f w ™* [N(E',9)-BG(E',9)] 
dn C J eff-ASl eff (E',0) ' 

(10) 

The individual ingredients will be discussed in detail in 
the following sections. 



target aluminum walls, negatively charged pions that are 
not separated from electrons by the PID cuts, and elec- 
trons originating from other processes, which are dom- 
inated by charge symmetric processes which produce 
equal numbers of positrons. Each of these potential back- 
grounds will be examined in the discussion that follows. 



1. Target Cell Backgrounds 

The quasielastic scattering from nucleons in aluminum 
nuclei can produce electrons at the same kinematics as 
those from elastic ep scattering. The scattering of the 
beam from front and back of the target cell wall produces 
backgrounds of this type which are difficult to isolate. 
Therefore, the corresponding background is determined 
by measuring the yield of events from a "dummy" tar- 
get, which is a mockup of the target ends. In order to 
minimize the data acquisition time, the total thickness 
of this dummy target was about 8 times the total cell 
wall thickness seen by the beam. After measuring the 
dummy yield, the total background from scattering from 
the target walls, BG W (E' ,9), was then determined from 



BG W (E', 



twQus 
tdQd 



N d (E',9)-C hT (E',9), 



(11) 



where Q w (d) is the total charge incident on the walls 
(dummy), t w id.) is the total thickness of the walls 
(dummy), and Nd{E' , 9) is the number of events collected 
for the dummy run after applying efficiency and deadtime 
corrections. 

The factor, Ch x (E', 9), corrects for the difference in ex- 
ternal bremsstrahlung emission due to the greater thick- 
ness of the dummy target. More precisely, this accounts 
for the fact that the distribution of events for thicker 
targets are more strongly shifted toward lower scattering 
energies (higher W) than those for thinner targets. The 
size of this correction was studied and was found to be 
less than a few tenths of a percent at all kinematics mea- 
sured, and typically less than 0.1%. Since this was the 
typical size of the uncertainty in this correction, we have 
taken Cb r = 1, and absorbed an additional 0.1%(0.1%) 
into the point-to-point (normalized) uncertainty in the 
background subtraction. 

The largest contribution to the uncertainties in the 
aluminum background subtraction comes from the uncer- 
tainties in the thickness of the cell wall of about 1.5% |l9j . 
However, the typical size of this background was on the 
order of 8% of the total yield, which leads to an uncer- 
tainty on the subtracted yield of only 0.12%. This un- 
certainty is approximately independent of the kinematics 
and run conditions. 



D. Backgrounds 

There are three physical processes that are possible 
sources of background counts to the elastically scattered 
electron yields. These are: electrons scattered from the 



2. Pion Backgrounds 

The rejection of negatively charged pions was accom- 
plished by placing requirements on both the number of 
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FIG. 9: Distribution of the number of photoelectrons collected 
in the Cerenkov detector for elastic kinematics of E — 2.2 GeV 
and p = 1 GeV/c. 




0.5 1 1.5 0.5 

Normalized Colorimeter Energy Deposit 



Cerenkov photoelectrons collected and the energy depo- 
sition of the particle in the calorimeter. The count distri- 
bution of photoelectrons collected in the HMS Cerenkov 
at an HMS momentum of 1 GeV is shown in Fig. [§] For 
electrons, this is a Poisson distribution, with a mean of 
approximately 10 photoelectrons. For pions, the number 
of photoelectrons produced should be zero. However, pi- 
ons can produce J-rays (electron knockout) in the ma- 
terials immediately preceding the Cerenkov detector and 
some of these "knock-on" electrons can produce Cerenkov 
radiation, with the probability of 5-ray production in- 
creasing with energy. With a requirement of more than 
2 photoelectrons, this decreases the pion rejection fac- 
tor from the maximum value of about 1000 : 1 found at 
low energies. However, this doesn't cause any significant 
pion contamination above this cut since the worst 7r~/e~ 
ratios are at low scattering energy where the rejection 
factor is the largest. 

The fractional energy deposition in the calorimeter, 
both before (unshaded region) and after (shaded region) 
applying the Cerenkov requirement of > 2 photoelectrons 
to select electrons, is shown in Fig. [H5| for the four kine- 
matics that exhibited the worst 7r~/e~ ratio. The frac- 
tional energy deposition of the particles is calculated by 
dividing the energy collected in a fiducial region about 
the track in the calorimeter by the momentum deter- 
mined from the track reconstruction. Even for these 
worst cases, it is evident that the Cerenkov requirement 
alone does a good job of removing pions. To further 
insure a clean electron sample, a requirement that the 
fractional energy deposited in the calorimeter be greater 
than 0.7 was also applied. The pion background after ap- 
plying both the Cerenkov and calorimeter requirements 
is estimated to be less than 0.1%. 

The Cerenkov efficiency, using a 2 photoelectron cut, 
was found to be 99.6%, independent of the energy. This 
is because the shape of the electron distribution does not 
depend on the particle's energy, resulting in the same 
fraction of electrons in the tail being removed. This is not 
true for the calorimeter, since a fixed energy resolution re- 
sults in an increase in the width of the electron fractional 
energy distribution at lower energies. The calorimeter 
cut efficiency decreases from a maximum of about 99.5% 
at energies above 3 GeV to about 98.5% at an energy of 
0.6 GeV. The run-to-run uncertainties on the efficiencies 
were estimated from Gaussian fits of the distribution of 
efficiencies determined for each run from the entire E94- 
110 elastic data set, and were found to be 0.1% for the 
Cerenkov detector and 0.1% for the calorimeter. 



E. Acceptance Corrections 



FIG. 10: Distribution of the fractional energy deposition in 
the calorimeter for events both before (open area) and af- 
ter (shaded area) applying a requirement of > 2 photoelec- 
trons in the Cerenkov detector. 



Whether a scattered electron reaches the detector stack 
or is stopped by hitting the edge of the collimator or one 
of the various apertures in the HMS magnet system and 
beam pipe is dependent upon several factors, including: 
1) the electron momentum, 2) the in-plane and out-of- 
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FIG. 11: HMS effective solid angle plotted in 8-5P/P space, 
using the binning described in the text. 

plane scattering angles, and 3) the vertex position. How- 
ever, the physics depends only upon the momentum and 
full scattering angle, 9 — cos -1 [cos(X') cos(6> c — Y')], so 
that for a fixed central spectrometer angle, 9 C , it is con- 
venient in what follows to consider only the E' and 9 
dependence of the acceptance averaged over the vertex 
coordinates. 

Using a model of the spectrometer, the fractional ac- 
ceptance, A(E',6), is calculated by generating Monte 
Carlo events and taking the ratio of the number of de- 
tected events to the number of generated events for each 
bin in phase space. That is 

A(E', 9) = N acc (E' gen ,9 gen )/N gen (E' gen ,9 gen ), (12) 

where N gen (E' gen , 9 gen ) is the number of events generated 
and N acc (E' gen ,9 gen ) is the number of events accepted 
in a given (E' gen ,9 gen ) bin. The gen subscripts denote 
that the kinematics used for the binning are as gener- 
ated. The fractional acceptance as defined here is simply 
a probability. However, it is evident that A(E',9) de- 
pends upon the solid angle, Afl gen (9) 1 into which events 
are generated. The "effective" solid angle coverage for 
each 2-dimensional bin is 

AQ eff {E',0) ee A(E',9)-An gen (A9), (13) 

and is independent of the size of Afl gen (A9). For exam- 
ple, increasing the generation limits of the out-of-plane 
angle, X' , from ±100 mrad to ±150 mrad will decrease 
A(E',0) since \X'\ > 100 mrad is already outside of the 
collimator aperture. However, Afl gen (A9) will increase 
accordingly and AQ e ff will remain unchanged. 



We note that the determination of A(E',9) does not 
require generating the events uniformly provided that the 
number generated in each part of phase space is known. 
The assumption here is that events generated in a given 
(E',6) bin are not detected in another (E',9) bin. The 
fractional acceptance as defined is then simply the prob- 
ability that an event generated in a given bin will be 
detected in that bin, and, therefore, the correction to 
the yield due to the fractional acceptance in each bin is 
1/A(E',9). If the bin-to-bin migration is small then it 
is already approximately accounted for by redefining the 
acceptance in Eauation ll2l to 

A{E', 9) = N acc (E' rec , 9 rec )/N gen (E' gen ,6 gen ), (14) 

where the rec subscripts denote the kinematics as re- 
constructed. The AQ e ff distribution extracted from the 
HMS model for E 1 = 2.8 GeV and 9 C = 12.5° is shown 
in figure ^] The shape in 9 — 9 C is dominated by the 
octagonal collimator, which largely determines the HMS 
solid angle acceptance. We note that the acceptance is 
not symmetric in the full scattering angle when the out- 
of-plane angle contributes significantly (i.e. at forward 
in-plane spectrometer angles) , even though the HMS has 
a high degree of symmetry about the in-plane scattering 
angle. This is because any out-of-plane angle will always 
result in a larger full scattering angle. 

The solid angle defined by the HMS collimator is about 
6.75 msr for a point target. This is slightly reduced for a 
4 cm extended target and the reduction becomes larger 
as the spectrometer is moved to larger angles. At the 
smallest angle measured of 9 C — 12.5°, the average solid 
angle acceptance due to the collimator for a momentum 
bite of \6p/p\ < 8% was determined from the HMS model 
to be 6.714 msr. The reduction due to all other apertures 
resulted in a further reduction of only 2.5% to 6.612 msr. 
At the largest angle measured of 9 C = 80°, the average 
solid angle acceptance due to the collimator was deter- 
mined to be 6.685 msr, with a further reduction due to 
other apertures of 5.2% to 6.335 msr. For this momentum 
bite the largest reduction of events after the collimator 
is in the second quadrupole. 

The normalization uncertainty on the acceptance cor- 
rections was estimated by combining in quadrature an 
uncertainty of 0.7% stemming from the reduction in solid 
angle due to apertures other than the collimator (more 
than one fourth the total at 9 C — 12.5°) and an uncer- 
tainty of 0.4% due to the modeling of the HMS optics. 

The optical properties of the HMS have been well stud- 
ied p6| utilizing several techniques and a large amount 
of dedicated optics data taken during many experiments 
over nearly a decade. For the HMS, the optical transport 
of charged particles through the spectrometer is indepen- 
dent of the momentum setting to a very high degree. The 
Afl distribution at a given 9 C is then only dependent on 
the energy setting through the dependence of the resolu- 
tion (including energy straggling) and multiple scattering 
effects in the spectrometer. 
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F. Elastic Peak Integration 



As already noted, the scattering energy distribution of 
the elastic peak at an individual 8 value is broadened 
from the 5 function expected in the OPEA due to sev- 
eral effects. These include energy resolution effects, and 
energy loss due to both ionization and bremsstrahlung 
emission. A typical peak distribution for a single 8 bin 
is shown in Fig. ^| The lower limit of integration is 
chosen to both minimize the loss of events due to res- 
olution smearing and to minimize the sensitivity to po- 
tential backgrounds, while the upper limit is chosen to 
include as much of the peak as possible yet to be below 
the threshold for inelastic it production at W 2 w 1.16 
GeV 2 . 

The sensitivities to both the lower and upper limits 
were studied and were found to be small. For the lower 
limit, the inscnsitivity indicates that the aluminum back- 
ground subtractions are correctly handled. For the up- 
per limit, it indicates that both the resolution effects 
and the shape of the bremsstrahlung distribution are ac- 
counted for reasonably well. Once the upper limit has 
been chosen, the fraction of the distribution that is out- 
side this limit is accounted for by the correction fac- 
tor RC(W 1 2 MJX ). If the effects of bremsstrahlung, energy 
straggling, and resolution are well understood, then a 
corresponding peak integration should be independent of 
the energy cutoff chosen, once the corresponding radia- 
tive correction has been applied. 

The sensitivity of the extracted cross section to the 
energy cutoff is shown in Fig. ^5] For each kinematic 
setting, the ratio of the integrated distribution for an en- 
ergy cutoff of W 2 iax to that for a cutoff of 1.15 GeV 2 is 
plotted versus e. In the upper plot W? nax — 1.10 GeV 2 , 
in the middle plot W^ ax = 1.05 GeV 2 , and for the bot- 
tom plot W^ax = 1-0 GeV 2 . The typical point-to-point 
difference in tail-corrected integration is less than 0.3% 
for the largest change in the cutoff value and shows little e 
dependence. We take this as the estimated random point- 
to-point uncertainty on this procedure and include it in 
the uncertainty of the radiative corrections. Addition- 
ally, we note that the normalization difference of about 
1% between between the smallest and largest values for 
Wmax i s hkely due to a combination of unoptimized res- 
olution matching and the approximate handling of the 
energy straggling in the simulation used for generating 
the acceptance corrections. However, this optimization 
becomes much less important as more of the peak is in- 
tegrated. We take 0.35% as the estimated normalization 
uncertainty on this procedure and include it in the un- 
certainty of the radiative corrections. 

Finally, we note that, except for the three measure- 
ments at beam energies below 2 GeV, the same W^ ax 
value was used for all kinematic settings. This was pos- 
sible because of the large acceptance of the HMS spec- 
trometer. This is in contrast to previous precision mea- 
surements d, , in which the spectrometer acceptance 
determined the maximum W 2 nrr at each kinematic set- 
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FIG. 12: Sample W 2 count distribution measured for elasti- 
cally scattered electrons in a single 9 bin after subtraction of 
the Al quasielastic contribution determined from the dummy 
target. 
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Bin-Centering and Averaging 



After performing the peak integration, the cross sec- 
tion is then extracted for each 8 bin. Often, the statistics 
taken in each bin are small (< 2000 counts). In order to 
improve the statistical accuracy, one would like to com- 
bine the data from all 8 bins. If the cross section did not 
depend (or depended only linearly) on the scattering an- 
gle, the cross sections extracted in each bin could simply 
be averaged. This is not the case, however. The HMS 
spectrometer has a relatively large acceptance of about 
±1.8 degrees in the scattering angle. Therefore, the cross 
section can vary greatly across the angular acceptance. 
At some kinematics, this variation can be a factor of 3 
or more (and strongly non-linear) across the acceptance. 
In order to average the cross sections in each 8 bin, the 
8 dependence of the cross section must be corrected for. 

This correction is called "8 bin-centering" (BC), and 
our prescription for it is straightforward. Since we would 
like to quote the cross section at the central angle of the 
spectrometer, the following correction is applied to each 
8 bin: 



da{8) 



dfl 



BC,i 



T Mod 



T Modt 



(15) 



where 6 is the central angle, 8i is the angle for the 
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FIG. 13: (Color online) Ratio of the cross section calculated 
with a cutoff in the tail integration of W 2 = 1.10 GeV 2 (top), 
W 2 = 1.05 GeV 2 (middle), and W 2 = 1.0 GeV 2 (bottom), 
relative to a cutoff of W 2 = 1.15 GeV 2 . 
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FIG. 14: Cross section extracted in each 6 bin across the HMS 
acceptance. The various data are discussed in the text. 



bin, and a is the value of a cross section model. For 
this procedure to be valid, care must be taken to sub- 
tract all backgrounds and to apply all corrections that 
have a 9 dependence, bin-by-bin. This includes radia- 
tive corrections. The bin-centered cross sections can then 
be averaged over the 9i to give the measured cross sec- 
tion at the central spectrometer angle. This was done 
as a weighted average, where the inverse of the square of 
the full statistical errors was used as a weighting factor. 
The statistical errors take into account the statistics of 
both the hydrogen and subtracted target endcap events, 
as well as the acceptance correction uncertainties due to 
statistical errors on the Monte Carlo generation. 

An example of this procedure is presented in Fig. [21 
Shown is the cross section extracted at a beam energy of 
3.12 GeV and a central HMS angle of 12.5°, before both 
acceptance and BC corrections (triangles). Also plot- 
ted is the cross section after applying acceptance cor- 
rections (squares) and after applying both acceptance 
and bin-centering corrections (circles). Only statistical 
uncertainties in the data are included in the error bars 
shown. However, the calculated acceptance corrections 
for bins at the edge of the acceptance can have large frac- 
tional errors, as they are very sensitive to both accurate 
modeling of the multiple scattering processes and small 
variations in the positions of apertures like the collima- 
tor. To minimize the effects of such sensitivities, bins at 



the edge of the 9 acceptance where the calculated accep- 
tance was below some minimum value were neglected in 
the averaging procedure. The angular acceptance limits 
used in the present experiment are represented by the 
vertical dashed lines in Fig. The cross sections ob- 
tained after averaging over these limits were found to 
be quite insensitive to the effects described above. Un- 
certainties associated with these effects were studied by 
adjusting aperture/target positions, magnet fields, and 
multiple scattering distributions within reasonable limits 
and determining the corresponding acceptance function. 
The cross section extracted with this acceptance function 
was typically found to agree with that using the nominal 
acceptance function to within 0.5%. 

The shape of the uncorrected distribution is the con- 
volution of the 9 dependence of the cross section with 
the acceptance of the spectrometer, with the latter being 
primarily defined by the collimator. Correcting for the 
acceptance leaves only the 9 dependence, which is then 
removed by the bin-centering corrections. The resulting, 
fully corrected, cross section should then be a constant 
across the 9 acceptance and equal to the cross section at 
the central 9 to within statistical fluctuations. This is 
indeed the case, allowing for small variations which are 
mostly due to an imperfect optics model for the spec- 
trometer, as well as the approximate treatment of the 
multiple scattering effects in the simulations. 
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(a) Vacuum polarization (b) Vertex corrections 





(c) Two photon exchange (d) Internal Bremsstrahlung 



FIG. 15: Feynman diagrams for higher order QED processes, 
including (a) vacuum polarization, (b) vertex corrections, (c) 
two photon exchange, and (d) internal bremsstrahlung emis- 
sion from electrons. Not shown are diagrams for proton 
bremsstrahlung, or external bremsstrahlung. 



H. Radiative Corrections 

While the form factors can be easily extracted from the 
Born cross section for single photon exchange, the cross 
section that is measured in a scattering experiment in- 
cludes higher order electromagnetic processes which are 
depicted in Fig. These processes can be categorized 
into two types: 1) the internal processes, which origi- 
nate due to the fields of the particles at the scattering 
vertex, and 2) the external processes, which originate 
due to the fields of particles in the bulk target materi- 
als. The internal processes include: vacuum polarization, 
vertex corrections, two-photon exchange, and (internal) 
bremsstrahlung emission in the field of the proton from 
which the scattering took place, while the external pro- 
cess is due to (external) bremsstrahlung in the field of a 
proton in the material either before or after the scattering 
vertex. 

The radiative correction factors, which account for 
these higher order processes, were calculated using the 
same procedure as the high precision SLAC data [1, 
and discussed in detail in Reference 0- T his p rocedure 
is based on the prescription of Mo and Tsai |2lll22T]. The 
radiative correction (RC) is applied to the measured cross 
section (after integration of the bremsstrahlung tail) as 
a multiplicative factor, i.e., 



da (6) 



RC- 



da{9) 

dn 



(16) 



Mi 



The corrections calculated include only the infrared di- 



vergent contributions from the two-photon exchange and 
proton vertex diagrams, with the non-divergent contri- 
butions having been previously estimated |22( to be less 
than 1%. However, the topic of two-photon exchange has 
recently been of renewed theoretical interest [30l l3ll l32] | 
in light of the discrepancy between the elastic form fac- 
tor ratios extracted from Rosenbluth separations of cross 
section measurements and those measured in polarization 
transfer experiments. In addition, a recent study |2^| of 
the world's data on the ratio of elastic cross sections for 
e + p to e~p has recently been made to look for evidence 
of two-photon effects. 

The radiative corrections applied at each kinematic 
setting are listed in Table ITTT1 We note that these correc- 
tion factors are significantly smaller than those applied in 
References |g, |9| . This is mostly due to integrating more 
of the radiative tail, but also to the reduction of external 
bremsstrahlung afforded by our much shorter target. 

The uncertainties in the radiative correction procedure 
were studied in References 0, @ and were estimated to 
be 0.5% point-to-point and 1.0% nomalized. To these we 
have added, in quadrature, the additional contributions 
resulting from the radiative tail integration, which are 
discussed in Section IVTH . 



I. Additional Corrections Applied to Large E' Data 

The four cross section measurements at scattering en- 
ergies greater than 3.5 GeV (indicated with a * in Ta- 
blc lIII|) required two additional corrections that were not 
needed for the rest of the data set. These corrections, 
which will be discussed in what follows, account for the 
following two effects: 1) an additional change in the ef- 
fective target density, and 2) a small misfocusing in the 
spectrometer optics, relative to nominal. 

During the data taking it was discovered that the fan 
controlling the flow of hydrogen through the cryogenic 
target had been inadvertently lowered from the 60 Hz 
nominal speed to 45 Hz. This resulted in larger local- 
ized boiling due to the beam and effectively lowered the 
density of hydrogen. To account for this, high statistics 
runs (« 0.1%) were performed at both fan speeds and 
the size of the effect was determined from the ratio of 
cross sections. Since the effect of target boiling was al- 
ready measured for a fan speed of 60 Hz, an additional 
correction to the yields of +0.6% ± 0.2% was included for 
the data taken at the lower speed. 

In addition to the correction for the difference in target 
density, the cross section measurements at large E' were 
also corrected for a slight misfocusing of the spectrome- 
ter. A TOSCA H| model of the HMS dipole indicated 
that it would start to suffer from saturation effects start- 
ing at currents corresponding to a momentum of 3.5 GeV 
for the central ray and that this effect would increase 
quadratically with momentum. This correction was in- 
cluded when setting the current in the dipole for this ex- 
periment, but was later found from both kinematic and 
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optics studies to not be needed. These studies indicated 
that the actual HMS field continues to increase linearly 
with current to a high accuracy up to the highest mo- 
mentum tested of 5.1 GeV. 

The missetting of the dipole field due to this un- 
needed saturation correction resulted in reconstructing 
the wrong scattering angle by an amount that varied 
across the angular acceptance, and led to a depletion 
of events that reconstructed in the region of the angular 
acceptance used in the analysis. The correction to the 
reconstructed scattering angle was determined by requir- 
ing W = M p for each 9 bin, and this correction was then 
fit as a function of 9 across the angular acceptance. The 
result of this fit was applied as a correction to the scatter- 
ing angle event-by-event, and effectively reshuffled events 
back into the depleted bins. The effect of this correction 
on the cross section was largest at the highest scattering 
momentum in this experiment of 4.7 GeV and resulted in 
a 1.6% increase. The uncertainty on this correction was 
estimated to be 0.4% and was assumed to be the same 
at all four kinematics where a correction was applied. 

VI. SYSTEMATIC UNCERTAINTIES 

The estimated systematic uncertainties for the exper- 
iment are listed in Table [I] for those that are assumed 
random point-to-point in e and in Table ITTI for those that 
effect the overall normalization uncertainty only. The 
quadrature sum of the point-to-point and normalization 
uncertainties gives the absolute uncertainties on the cross 
section measurements. The point-to-point uncertainties 
are those that depend upon variable run conditions or 
kinematics. Discussions of the uncertainties presented 
here can be found in earlier sections of the text. 

TABLE I: E94-110 point-to-point systematic uncertainties. 



Experimental Quantity Uncertainty Aa/a 



(pt-pt) 


Beam Energy 


4 ■ 1(T 4 


0.0024 


Scattering Angle 


0.2 mrad 


0.0026 


Target Density 


0.1% 


0.001 


Target Length 


0.1% 


0.001 


Beam Charge 


0.2% 


0.002 


Acceptance 


0.5% 


0.005 


Detector Efficiency 


0.15% 


0.0015 


Tracking Efficiency 


0.25% 


0.0025 


Deadtime Corrections 


0.14% 


0.0014 


Target Cell Background 


0.2% 


0.002 


Radiative Corrections 


0.6% 


0.006 




Total 


0.0097 



VII. RESULTS 



TABLE II: E94-110 normalization uncertainties. 



Experimental Quantity Uncertainty Aa/a 



(Norm) 


Beam Energy 


4 • 10" 4 


0.0024 


Scattering Angle 


0.4 mrad 


0.0053 


Target Density 


0.4% 


0.004 


Target Length 


0.3% 


0.003 


Beam Charge 


0.33% 


0.0033 


Acceptance 


0.8% 


0.008 


Detector Efficiency 


0.4% 


0.004 


Tracking Efficiency 


0.3% 


0.003 


Deadtime Corrections 


0.1% 


0.001 


Target Cell Background 


0.3% 


0.003 


Radiative Corrections 


1.1% 


0.011 




Total 


0.017 



The complete set of Born cross sections extracted from 
these precision elastic scattering measurements are listed 
in Appendix A. In total, measurements at 28 different 
kinematics were included in the final data set, covering 
a range in Q 2 from approximately 0.4 GeV to 5.5 GeV. 
In addition, a significant range in e was covered, even 
though the elastic kinematics were not specifically opti- 
mized for Rosenbluth separations. 



A. Comparisons to Fits of the Existing World's 
Data 

Comparisons of these cross sections were made to re- 
cent fits of the previous world's data set. These included 
the fit of Brash et al. from and the fits of Arrington 
from [l2^ . Both analyzes performed a fit to the com- 
bined cross section and polarization transfer data, while 
that of Arrington also included a fit to the cross section 
data alone. It should be stressed that, although nearly 
the same data sets were included in the fits, the method 
of combining the polarization transfer data differed. 

In the work of Brash et al., the polarization transfer 
results for the Q 2 dependence of Ge/Gm were used as a 
constraint to do a one-parameter refit of the Rosenbluth 
data, from which Gm was extracted. A \ 2 minimization 
was then performed to fit Gm as a function of Q 2 . In 
the work of Arrington, two fits were performed: one that 
included the polarization transfer data, and one that did 
not. However, here the polarization transfer results for 
Ge/Gm were included with the cross section measure- 
ments as equally weighted points in a x 2 minimization to 
fit both form factors simultaneously. This follows closely 
the older work of Walker Q on global fits to cross section 
data. 

The ratio of the data to the Arrington fit of cross sec- 
tion data is shown in Fig. ^] versus 9, e, and E' . The 
inner error bars represent the purely statistical uncertain- 
ties, while the full error bars include the point-to-point 
uncertainties as well. This fit is observed to describe the 
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TABLE III: Table of measured ep elastic cross sections. The systematic uncertainties listed are the estimated point-to-point 
uncertainties for each kinematic setting. In addition there is scale uncertainty of 1.7%. 





e 




Q 2 




e 


a 


Act (stat) 


Act (sys) 


RC 


(GeV) 


(D 


eg) 


(GeV/c) 2 






(ph/sr) 


(/ttb/sr) 


Ob/ sr ) 






1.148 


47 


.97 


0. 


.6200 





.6824 


0.1734E-01 


0.35E-04 


0.16E-03 


1. 


,055 


1.148 


59 


.99 


0. 


.8172 





.5492 


0.4813E-02 


0.11E-04 


0.45E-04 


1 


,047 


1.882 


33 


.95 


0. 


.8995 





.8104 


0.1464E-01 


0.28E-04 


0.14E-03 


1 


,098 


2.235 


21 


.97 


o 


,6182 





.9187 


0.1098E+00 


0.40E-03 


0.11E-02 


1 


,120 


2.235 


31 


.95 


1 


1117 





.8226 


0.8938E-02 


0.25E-04 


0.86E-04 


1 


,121 


2.235 


42 


.97 


1 


,6348 





.6879 


0.1184E-02 


0.54E-05 


0.11E-04 


1 


,114 


2.235 


58 


.97 


2 


2466 





.4885 


0.1522E-03 


0.86E-06 


0.14E-05 


1 


,100 


2.235 


79 


.97 


2 


,7802 





.2843 


0.2868E-04 


0.21E-06 


0.27E-06 


1 


,083 


3.114 


12 


.47 





,4241 





.9740 


0.8968E+00 


0.25E-02 


0.96E-02 


1 


,150 


3.114 


15 


.97 





,6633 





.9553 


0.1892E+00 


0.80E-03 


0.20E-02 


1 


,151 


3.114 


19 


.46 





,9312 





.9308 


0.4802E-01 


0.12E-03 


0.49E-03 


1 


,154 


3.114 


32 


.97 


2. 


.0354 





.7835 


0.1017E-02 


0.26E-05 


0.99E-05 


1 


,154 


3.114 


40 


.97 


2 


6205 





.6726 


0.2125E-03 


0.74E-06 


0.20E-05 


1 


,150 


3.114 


49 


.97 


3 


.1685 





.5480 


0.5568E-04 


0.28E-06 


0.53E-06 


1 


,147 


3.114 


61 


.97 


3 


,7261 





.4026 


0.1487E-04 


0.10E-06 


0.14E-06 


1 


,137 


3.114 


77 


.97 


4. 


,2330 





.2574 


0.4260E-05 


0.49E-07 


0.40E-07 


1 


,121 


4.104 


38 


.97 


3 


,7981 





.6578 


0.4919E-04 


0.31E-06 


0.47E-06 


1 


,185 


4.104 


45 


.96 


4 


,4004 





.5528 


0.1582E-04 


0.29E-06 


0.15E-06 


1 


,181 


4.413 


44 


.98 


4. 


,7957 





.5526 


0.1098E-04 


0.20E-06 


0.11E-06 


1 


,192 


4.413 


50 


.99 


5 


,2612 





.4686 


0.4898E-05 


0.73E-07 


0.47E-07 


1 


,188 


★5.494 


12 


.99 


1. 


,3428 





.9655 


0.4055E-01 


0.11E-03 


0.50E-03 


1 


,214 


★5.494 


17 


.96 


2 


,2878 





.9239 


0.2866E-02 


0.13E-04 


0.33E-04 


1 


,220 


★5.494 


20 


.47 


2 


,7822 





.8955 


0.9824E-03 


0.41E-05 


0.11E-04 


1 


,223 


★5.494 


22 


.97 


3 


.2682 





.8627 


0.3770E-03 


0.15E-05 


0.42E-05 


1 


,226 


5.494 


25 


.47 


3 


,7385 





.8261 


0.1608E-03 


0.14E-05 


0.18E-05 


1 


,226 


5.494 


27 


.97 


4 


,1867 





.7865 


0.7749E-04 


0.61E-06 


0.77E-06 


1 


,228 


5.494 


32 


.97 


5, 


,0031 





.7023 


0.2149E-04 


0.23E-06 


0.21E-06 


1 


,226 


5.494 


35 


.48 


5 


,3699 





.6593 


0.1267E-04 


0.20E-06 


0.12E-06 


1 


,225 



data set very well over the entire kinematic range. The 
X 2 (calculated using only the point-to-point uncertain- 
ties and after removing the average 0.6% normalization 
difference between the current data set and the previ- 
ous world data set of cross section measurements) dis- 
tribution was found to be well described by a Gaussian 
distribution with a width corresponding to an average 
uncertainty of about 1.0%, consistent with the estimated 
errors combining the systematic point-to-point and sta- 
tistical uncertainties in quadrature. For each of the three 
fits previously described, the total \ 2 P er degree of free- 
dom (xl) to the data was calculated. The results for the 
region above Q 2 = 1 (GcV/c) 2 , where the discrepancy 
between the cross section and polarization results differ 
significantly, were found to be xt = 0-76 (Arrington fit to 
cross sections), 1.06 (Arrington fit including polarization 
transfer results), and 2.95 (Brash et al fit), allowing the 
overall normalization to vary. 

These results are interesting for two reasons. Firstly, 
the full data set above Q 2 = 1 (GcV/c) 2 favor the fit to 
cross section data only over the fits that includes the po- 
larization transfer data. Secondly, the data favor the Ar- 
rington prescription for combining the cross section and 
polarization transfer data over the prescription of Brash 
et al. This does not resolve the inconsistency between the 



Rosenbluth and polarization transfer results, but rather 
underscores a consistency in the global cross section data 
set including these new measurements. The discrepancy 
with these and the polarization transfer measurements is 
highlighted further. 

B. Rosenbluth Extractions of Form Factors 

The individual Sachs form factors were extracted from 
the cross section data at seven different Q 2 values via 
the Rosenbluth separation method. This required that 
the cross section measurements at similar Q 2 values be 
grouped. Since none of the measurements were taken at 
precisely the same Q 2 , a correction factor was applied to 
some of the cross sections in each group to evolve to a 
common Q 2 . The correction factor for this Q 2 evolution 
was calculated from fits to previous data via 

da(Q 2 c ,e) da(Q 2 ,e) <r M ° d (Q 2 c ,e) 

dn dfl a M ° d {Q 2 ,e)' { ' 

where <j Mod is the value given by the fit of Arrington to 
the cross sections data, and Q 2 , Q 2 , represent the values 
before and after the evolution, respectively. 
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FIG. 16: (Color online) The ratio of measured ep elastic cross 
sections to the Arrington fit of previous cross section data, 
as a function of scattering angle, longitudinal photon polar- 
ization, and scattering energy. The solid line indicates the 
average ratio of 1.006 and the dashed lines indicate ±2% of 
this value. 
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FIG. 17: Rosenbluth separations of the form factors. Plotted 
is the reduced cross sections (xlOOO) versus e. 



In order to perform Rosenbluth separations at a par- 
ticular Q 2 , the following two conditions on the data were 
required: 1) each separation must contain three distinct 
e points, and 2) the Q 2 evolution for each s point must 
constitute less than a 15% correction. The sensitivity 
of the extracted form factors on the model used for the 
Q 2 evolution was found to be much less than the uncer- 
tainties. Plots of the reduced cross sections versus e are 
presented in Fig.E|f or each Q 2 . Also presented are the 
results of the linear fit. The error bars on each point 
represent the total point-to-point uncertainties, includ- 
ing both statistical and systematic uncertainties added 
in quadrature. 



The results for the ratio /iGe/Gm extracted from the 
current data set are presented in Fig. 1181 along with pre- 
vious extractions from both cross section [8j and polar- 
ization transfer 0, Q data. The current data are seen 
to agree well with the previous cross section data, while 
being in significant disagreement with the polarization 
transfer results. The error bars on each point represent 
the uncertainties obtained for the fit parameters, while 
the hatched band at the top of the figure represents that 
due to the estimated 0.4 mrad uncertainty in the abso- 
lute scattering angle. To a large degree, an error in the 
scattering angle would shift the entire set of ratios up 
or down but would not significantly alter the trend ver- 
sus Q 2 , which is significantly different from that of the 
polarization transfer data. 
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TABLE IV: Table of the Rosenbluth extracted Sachs form factors relative to the dipole form factor, Gd = n^Q^m n)^ (Q 2 m 
(GeV/c) 2 . 



Q 2 

(GeV/c) 2 


GMp/ifJ-pGdip) 


Ge p / Gdip 


HGEp/GMp 


0.65 


0.968 ± 0.032 


1.035 ± 0.052 


1.069 ± 0.085 


0.91 


1.028 ± 0.019 


0.954 ± 0.053 


0.928 ± 0.067 


2.20 


1.050 ± 0.016 


0.923 ± 0.121 


0.878 ± 0.125 


2.75 


1.055 ± 0.010 


0.888 ± 0.114 


0.841 ± 0.109 


3.75 


1.044 ± 0.015 


0.873 ± 0.232 


0.837 ± 0.220 


4.20 


1.012 ± 0.012 


1.255 ± 0.157 


1.240 ± 0.163 


5.20 


1.007 ± 0.032 


1.183 ± 0.511 


1.176 ± 0.552 




0.2 -k - Walker Global 

- ■ - Polarization Transfer 
- • - This Experiment 
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1 



2 



Q 2 [(GeV/c) 2 ] 



between Rosenbluth and polarization transfer extractions 
of the ratio ^lGe/Gm- This data set shows good agree- 
ment with previous cross section measurements, indicat- 
ing that if a here-to-fore unknown systematic error does 
exist in the cross section measurements then it is intrinsic 
to all such measurements. 

A likely candidate, which has received much theoreti- 
cal interest recently |3(J, l3ll l32| , is possible contributions 
from two-photon exchange, which are not fully accounted 
for in the standard radiative corrections procedure of 
Mo-Tsai. Although it is currently unclear whether such 
an effect can fully explain the discrepancy, considerable 
progress is being made. 

Complementary to this theoretical effort is the recently 
completed experiment |33j in JLab Hall A which utilizes 
the so-called 'Super-Rosenbluth' technique to extract the 
form factor ratio. This experiment measured the pro- 
ton cross sections and is therefore sensitive to a different 
set of systematic uncertainties than the previous elec- 
tron cross section data. However, these measurements 
are still as sensitive to two-photon exchange effects as 
electron cross section measurements and will, therefore, 
provide a vital clue whether such effects are present. In 
any event, it is critical that the source of the discrepancy 
be found if there is to be any hope of extracting the Q 2 
dependence of the individual form factors. 



FIG. 18: (Color online) Extracted values of [iGe/Gm. Also 
shown are those extracted from previous Rosenbluth and po- 
larization transfer measurements. The error bars on the data 
include the full point-to-point uncertainties, while the hatched 
band at the top indicates the uncertainty due to the absolute 
uncertainty in the scattering angle of 0.4 mrad. 

VIII. CONCLUSION 

We have performed high precision measurements of the 
ep elastic cross section covering a considerable amount of 
the Q 2 — e space for which there exists a large discrepancy 
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